Freeway Ramp Metering Control 
Made Easy and Efficient 
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Abstract: "Model- free" control and the related "intelligent" proportional-integral (PI) con- 
trollers are successfully applied to freeway ramp metering control. Implementing moreover the 
corresponding control strategy is straightforward. Numerical simulations on the other hand need 
the identification of quite complex quantities like the free flow speed and the critical density. 
This is achieved thanks to new estimation techniques where the differentiation of noisy signals 
plays a key role. Several excellent computer simulations are provided and analyzed. 

Keywords: Traffic control, ramp metering, congestion, fundamental diagram, free flow speed, 
critical density, model-free control, intelligent PI controllers, identification, estimation, 
numerical differentiation. 
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(1) the control of freeway ramp metering, 

(2) the estimation of the free-flow speed and of the critical Computer experiments show that our control strategy 
density behaves better than ALINEA, which was a most remark- 
able breakthrough when introduced more than twenty 

Freeway ramp metering control, which should alleviate yea rs ago (Hadj-Salem, Blosseville, Davee & Papageorgiou 

congestions, is achieved via model-free control (Fliess & (1988); Hadj-Salem, Blosseville & Papageorgiou (1990); 

Join (2008, 2009)). It yields an intelligent proportional- Papageorgiou, Hadj-Salem & Blosseville (1991)). 2 Despite 

integral, or iPI, controller which a huge academic literature, which utilizes most of the 

• regulates the traffic flow in a most efficient way, cxistin S methods of modern control theory, whether with 
. is robust with respect to quite strong disturbances, lum P ed or wlth distributed parameter systems, ALINEA, 
. is easy to tunc and to implement, wmch 18 exploited in France and in many other countries, 

• does not need any precise mathematical modeling. 

ALINEA is an acronym of Asservissement LIN eaire d'Entree 

Remark 1.1. Model- free control, although quite new, has Autoroutiere. 

already been successfully employed in many concrete sit- ~See, e.g., (Papageorgiou, Diakaki, Dinopoulou, Kostialos & Wang 
uations: (2003); Smaragdis & Papageorgiou (2004)) for recent developments. 



remains until today to the best of our knowledge the 
only feedback-control law for ramp metering that has been 
implemented in practice. 

Computer simulations, on the other hand, need some 
kind of precise mathematical macroscopic modeling. They 
become therefore more subtle and complex than model- 
free control. This necessity yields a dichotomy which is 
analyzed in this paper for the first time. We utilize here 
ordinary differential equations, i.e., a macroscopic model 
of order two, due to Payne (1971) and improved by Pa- 
pagcorgiou, Blossevillc & Hadj-Salcm (1990). The corre- 
sponding model properties are quite sensitive to parameter 
variations and uncertainties. The free-flow speed and the 
critical density are estimated here via the fundamental 
diagram due to May (1990) thanks to recent differentiation 
techniques of noisy signals (Fliess, Join & Sira-Ramfrcz 
(2008); Mboup, Join & Fliess (2009)). Most of the existing 
methods for achieving real-time estimation employ in one 
way or the other the Kalman filtering (see, e.g., Mihaylova, 
Boel & Hegyi (2009); Wang & Papageorgiou (2005); Wang, 
Papageorgiou & Messmer (2008)). Their computational 
burden seems however quite higher than ours. 

Our paper is organized as follows. Model-free control and 
intelligent PI controllers are presented in Section 2. 3 Sec- 
tion 3 studies the application to a concrete example of an 
isolated ramp metering. After reviewing the identification 
techniques which arc connected to the fundamental dia- 
gram, important parameters corresponding to the same 
freeway are estimated in Section 4. Convincing computer 
simulations are also analyzed in Sections 3 and 4. Some 
concluding remarks are discussed in Section 5. 

2. MODEL-FREE CONTROL: A SHORT REVIEW 
2.1 Basics 

We restrict ourselves for simplicity's sake to a SISO system 
6, with a single input u and a single output y. We do 
not know any global mathematical description of &. We 
replace it by a "phenomenological" model, which is 

• valid during a short time lapse, 

• said to be ultra-local, 



,0) 



F 



(1) 



where 



the differentiation order v of y, which is 

■ chosen by the practitioner, 

• generally equal to 1, 
has no connection with the unknown differentiation 
order of y in ©; 

the constant parameter a has no a priori precise 

numerical value. It is determined by the practitioner 

in such a way that the numerical values of era and 

j/"' are of equivalent magnitude; 

F, which contains all the "structural" information, 

depends on all the system variables including the 

perturbations. 



2.2 Intelligent PI controllers 

Assume that we have a "good" estimate 4 [F) e of F and, 
for simplicity's sake, that v = 1 in Equation (1). 5 The de- 
sired behavior is obtained via an intelligent proportional- 
integral, or iPI, controller 




(2) 



where 



• y* is the output reference trajectory, 

• e = y — y* is the tracking error, 

• Kp, Kj are the usual gains. 

If Kj = 0, we have an intelligent proportional, or iP, 
controller: 




(3) 



Remark 2.1. Contrary to the situations with classic PI 
controllers, controllers (2) and (3) are easy to tune: they 
stabilize a pure integrator. 

Remark 2.2. See d'Andrea-Novel, Fliess, Join, Mounier & 
Steux (2010) for the explanation of the strange ubiquity 
of classic PIDs via the above viewpoint. 6 

2.3 Estimation of F 

Estimation of y If v = 1 in Equation (1), [F] e may be 
obtained via the estimate of y. Elementary differentiation 
filters do suffice in this situation where the sampling is 
rather crude. 

Another technique Rewrite Equation (2) as 
F = —au + y* — Kpe — Ki I e 



Corrupting noises are attenuated by integrating both sides 
on a short time interval. 7 It yields: 



-^approx ^ I 



i J (-au + ij* - K P e - A"/ J e^j dr (4) 

where -Fkpprox is a piecewise constant approximation of 
F. Equation (4) may be easily implemented as a discrete 
linear filter. 



3. FREEWAY RAMP METERING PRINCIPLE 

3.1 Generalities 

Consider the simple example of the freeway section de- 
picted in Fig. 1: 



• q r , in veh/h, is the ramp flow related to the control 
variable r <E {r min , r max }, 8 by q r = rq r , where 



mm 



-,Q sa t min r, 



Ps-Pc 



is the flow; 



See Fliess, Join & Riachy (2011) for a complete presentation. 



4 See Section 2.3. 

5 v = 1 is an appropriate choice for most of the concrete examples. 
See Fliess, Join & Riachy (2011) for an explanation. 

6 See also Fliess, Join & Riachy (2011). 

7 See Fliess (2006) for a mathematical explanation. 

8 Just as in Hegyi, Dc Schutter & Hellendoor (2005), we set r = 1 
for an unmetered on-ramp. 
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Figure 1. Freeway ramp metering principle 

• w represents the queue length in vehicles, 

• Qsat, is the on-ramp capacity in veh/h, 

• pmaxi Pc are respectively the maximum and the 
critical density. 

The ramp metering, or admissible control, consists to act 
on the traffic demand at the on-ramp origin in order to 
maintain the traffic flow in the mainstream section close 
to the critical density. 9 

Ramp metering strategies may be local (isolated ramp 
metering) or coordinated (Smaragdis & Papageorgiou 
(2004)). Isolated ramp metering makes use of real-time 
traffic measurements in the vicinity of each controlled 
on-ramp in order to calculate the corresponding suitable 
ramp metering flows. Coordinated ramp metering exploits 
the all available measurements of the considered portions 
of controlled freeway. We focus here on isolated ramp 
metering. 

3.2 Model- free ramp metering 

For the studied freeway section (see Fig. 1), Equation (1) 
becomes 10 



p 8 {t) = F(t) - ar{t) 



(5) 



The control variable r(t) is given via the intelligent con- 
troller iPI (2): 



r(t) 



1 

a 



-in 



K P e + Ki 



(0) 



where 



- p* is the reference trajectory. 

- e = p s — p* is the tracking error. 

The estimation of F is provided thanks to the following 
expression: 



[F(k)] e = [p s (k)] e - ar(k - 1)} 



where 

- k is the sampled time, 

9 The traffic demand is assumed to be independent of any control 
actions (see, e.g., Papageorgiou, Blosseville & Hadj-Salem (1990); 
Kostialos, Papageorgiou & Middelham (2001)). 

10 The traffic occupancy measurements are utilized for practical 
implementation. 
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Figure 2. General scheme of the input/output control 
- [»] e indicates an estimate of •. 11 

3.3 Implementation issues 

Implementing our model-free control and the related iPI 
controller is straightforward: 12 

• The gains Kp and Kj are easily tuned thanks to the 
first order ultra-local model (5). 

• The remarkable robustness properties follow from the 
excellent estimation [F] e of F. 13 

• the generation of the desired trajectory (density) p* 
is achieved thanks to the following algorithm: 

• Let Vfiiterod be the filtered mean speed and 
^threshold the speed threshold. 14 

• PdOi Pine, Pdoc denote respectively the initial den- 
sity, the increment and decrement of the desired 
density. 

■ If ^filtered > ^threshold, then p* = p t \ + pi nc . 

■ If ^filtered < ^threshold, then p* = PdO — Pdec- 

3- 4 Simulation results 

Our computer simulations are based on numerical data 
which are collected from the French freeway A4Y with one 
on-ramp (see Fig. 3 and Fig. 4). The software METANET 
(Papageorgiou (1983)) is utilized. 15 Although the mea- 
surements, i.e., the traffic volume in veh/h, are quite poor 
and noisy, the performances of our iPI controller (Fig. 7) 
are good. Congestions are alleviated as soon as they appear 
(Fig. 5 and Fig. 6). 

4. TRAFFIC FLOW PARAMETRIC ESTIMATION 

4- 1 Generalities 

The macroscopic models which are used for simulation 
purposes, are not only heuristic but also quite sensitive to 
parameter variations and uncertainties. The only available 
accurate physical law is the conservation equation. All 
other equations (speed equation and the fundamental di- 
agrams, for instance), are based on empirical observations 

11 See Section 2.3. 

12 See Fig. 2 for a corresponding block diagram scheme representa- 
tion. 

13 See Section 2.3. 

14 Concrete studies (see, e.g., Cete Mediterranee - Les etudes (2006)) 
have demonstrated that the level of service is highly degraded and 
the congestion phenomenon is at its maximum, when the mean speed 
of individual vehicle is about 30 km/h. The threshold of discomfort 
is reached, when this speed is equal to 85 km/h (see e.g. Cete 
Mediterranee - Les etudes (2006)). 
15 It is based on a second order macroscopic model. 



Figure 3. Aerial picture of the studied site (Source DiRIF) 
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Figure 4. Traffic demands: (-) mainstream, (-) on-ramp 



Figure 6. Mean speeds evolutions: (-) no-control case, (- 
control case 



100 




Figure 5. Traffic densities evolutions: (-) no-control case, 
(-) control case, (-) p* 

and coarse approximations. The main parameters such as 
the critical density and the free-flow speed are moreover 
subject to variations. 

4-2 Fundamental diagram 

The fundamental diagram due to May (1990) is given by 




Figure 7. Control variable 




where 



• pi is the density of the segment i, 

• V is the corresponding the mean speed, 



(7) 




Figure 8. Fundamental diagram example 

• Vf is the free-flow speed, 

• p c the critical density, 

• a is a model parameter. 

Although it is of heuristic nature (see Fig. 8), i.e., not 
derived from physical laws, it provides important param- 
eters for the macroscopic modeling which is used for the 
numerical simulations and to identify congestion and fluid 
zones. 

4-3 Identification 

The existing results in signal processing (see Fliess, 
Mboup, Mounier & Sira- Ramirez (2003); Fliess (2008); 
Mboup (2009)) should be extended to the arbitrary ex- 
ponent a in Equation (7). 

New setting Rewrite Equation (7) in the following form 
V( Pi ) = v f ex P [-Kp a ] (8) 
where K = The equality 



V 

Pi 



(10) 



The numerical derivation of noisy signals, developed in 
(Fliess, Join & Sira- Ramirez (2008); Mboup, Join & Fliess 
(2009)), has been already successfully implemented in 
many concrete applications (see, for example in intelligent 
transportation systems, Menhour, d'Andrea-Novel, Bous- 
sard, Fliess & Mounier (2011); Menhour, d'Andrea-Novel, 
Fliess & Mounier (2012); Villagra, d'Andrea-Novel, Fliess 
& Mounier (2009, 2011)). In order to summarize the gen- 
eral principles, let us start with the first degree polynomial 
time function pi(t) — oq + ait, t > 0, <Xo,ai G K. Rewrite 
it thanks to classic operational calculus (see, e.g., Yosida 
(1984)) pi as Pi = 2a + % Multiply both sides by s 2 : 



s Pi = aos + ai 



(11) 



Take the derivative of both sides with respect to s, which 
corresponds in the time domain to the multiplication by 
-t: 



s 2 —± + 2sP 1 
as 



ao 



(12) 



The coefficients oo,ai are obtained via the triangular 
system of equations (11)-(12). We get rid of the time 
derivatives, i.e., of sP]_, s 2 P\, and s 2 ^j^-, by multiplying 
both sides of Equations (11)-(12) by s~ n , n > 2. The 
corresponding iterated time integrals are low pass filters 
which attenuate the corrupting noises, which are viewed 
as highly fluctuating phenomena (Fliess (2006)). A quite 
short time window is sufficient for obtaining accurate 
values of ao, a±. 

The extension to polynomial functions of higher degree 
is straightforward. For derivatives estimates up to some 
finite order of a given smooth function / : [0, +oo) — y R, 
take a suitable truncated Taylor expansion around a given 
time instant to, and apply the previous computations. 
Resetting and utilizing sliding time windows permit to 
estimate derivatives of various orders at any sampled time 
instant. 



shows that p c may be deduced at once from a and K. 

Notation: If G is a function of pi, write G Pi its derivative 
with respect to pi. 

Write W the logarithmic derivative of V with respect to 
Pi- 

Vn 

W=-^ = -Kap a ~ x 

Thus 



(9) 



a- 1 



W Pl 

The identifiability of a follows at once. Equation (9) 
provides K and Equation (8) Vf. 

Remark £.1. Note that the second order derivative V p 2 of 

V is needed. 

Derivation with respect to time Consider pi and, then, 

V as functions of time t. The time derivatives are obtained 
using the following expression: 



4-4 Computer experiments 

The measurements concern the evolution of the mean 
speed and of the traffic occupancy depicted in Fig. 9. 16 
The data are 

• provided during five days with a sampling period of 
20 seconds, 

• quite poor and noisy. 

The modeling approximation, 17 and the numerical singu- 
larities, which are unavoidable in such a real-time setting, 
explain why our estimates do fluctuate to some extent. The 
results depicted in Fig. 10 do however show a satisfactory 
"practical" convergence towards values which are suitable 
for our simulation purposes. 



16 The occupancy measurements arc transformed into traffic density 
for simulations purposes. 

17 As already stated in Section 4.2 the fundamental diagram is only 
heuristic. 
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(a) Density evolution p vs time in seconds 



(a) Critical density estimation p c vs time in seconds 
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(b) Speed evolution v vs time in seconds 

Figure 9. Measured variables 

5. CONCLUSION 

The control design described in this paper will soon be 
implemented in practice. The pending patent 18 prevents 
us unfortunately from discussing future developments in 
traffic control. 

From a theoretical standpoint two major points might 
however be stressed: 

(1) The control of freeway ramp metering and of hydro- 
electric power plants is approached almost exclusively 
in the existing academic literature via a rather com- 
plex modeling where partial differential equations are 
often utilized. 19 The present work confirms what has 
already been obtained for hydroelectric power plants 
by Join, Robert & Fliess (2010), namely that an 




(b) Free-flow speed estimation v vs time in seconds 




18 See the acknowledgement below. 

19 See, e.g., some references in Aboua'issa, Fliess, Iordanova & Join 
(2011b). 



(c) Estimation of a vs time in seconds 

Figure 10. Estimated parameters 



elementary model-free control is enough for obtaining 
excellent results. 20 
(2) Computer simulations do necessitate a quite realistic 
modeling which implies a more subtle mathematical 
setting than model-free control. This fact which is 
stressed here for the first time will be further studied 
in the future. It might lead to a profound epistemo- 
logical revolution in applied mathematics and, more 
generally, in applied sciences, the consequences of 
which are not yet clear. 
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